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The  purpose  of  this  work  is  to  formulate  a general  theory  for  the 
dynamic  response  of  linear  damped  continuous  structural  members  subjected 
to  arbitrary  excitation  forces.  Typical  members  which  are  included  in 
this  theory  are  rods,  beams,  plates,  shells,  and  rotating  shafts.  The 
damping  models  to  be  considered  include  the  traditional  viscous  dashpot, 
so-called  hysteretic  (or  material  or  structural)  damping,  and  viscoelastic 
material  damping. 

A general  technique  for  finding  the  dynamic  response  is  modal 
analysis,  and  it  will  be  used  exclusively  in  the  following  developments. 
This  method  has  been  employed  extensively  for  both  undamped  discrete  and 


undamped  continuous  linear  systems  (see  for  example,  Hurty  and  Rubinstein 
[1]^"  and  Meirovitch  [2]).  However,  when  damping  is  present,  the  use  of 
modal  analysis  has  been  largely  restricted  to  self-adjoint  systems  of 
equations  which  may  be  uncoupled  with  the  classical  or  undamped  normal 
modes.  Meirovitch  [2]  presents  several  examples  of  such  systems  which 
include  structural  members  with  proportional  viscous  and  linear  hysteretic 
damping.  For  viscoelastic  structural  members,  Valanis  [3]  was  the  first 
to  propose  a suitable  modal  solution.  By  assuming  that  Poisson's  ratio 
remains  constant,  the  dynamic  viscoelastic  problem  may  be  resolved  into 
a quasi-static  viscoelastic  problem  and  a dynamic  undamped  elastic 
problem.  Robertson  and  Thomas  [4]  used  a similar  technique  which  allowed 
more  general  boundary  and  non-zero  initial  conditions  to  be  considered. 


lumbers  in  bracket 


s refer  to  References  at  end  of  paper. 
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an  example  of  this  technique  as  applied  to  a self-adjoint  viscoelastic 
Timoshenko  beam  is  given  by  Robertson  [5]. 

When  other  forms  of  damping,  such  as  non-proportional  viscous 
damping,  are  included  in  a structural  member  model  the  classical  modes 
cannot  uncouple  the  equations  of  motion.  Hence,  modal  analysis  in  the 

"s. 

classical  sense  may  not  be  used  to  determine  the  response.  However,  for 
a discrete  system  with  viscous  damping,  it  has  been  shown  by  Foss  [6] 
that  a resolution  of  the  response  is  possible  by  using  a new  set  of 
orthogonal  eigenfunctions.  These  are  the  so-called  damped  modes,  which 
can  be  shown  to  satisfy  a different  form  of  orthogonality  [6].  In  a 
similar  manner,  the  response  of  a continuous  member  with  non-proportional 
viscous  damping  was  determined  by  O' Kelley  [8]  by  using  the  damped  mode 
shapes.  In  addition,  Caughey  and  O’Kelley  [8]  studied  the  limitations 
of  modal  analysis  for  this  kind  of  general  viscous  damping  by  determining 
necessary  and  sufficient  conditions  for  which  the  undamped  modes  can  be 
used  to  uncouple  the  dynamic  motion  (i.e.  proportional  viscous  damping). 
However,  few  applications  have  appeared  in  the  literature  applying 
damped  modal  analysis  to  continuous  members  with  viscous  damping.  For 
example,  Lund  [9]  used  damped  eigenfunctions  to  find  the  dynamic  response 
of  a rotating  shaft  on  damped  bearings.  The  approach  taken  was  not 
general  and  it  involved  a great  deal  of  tedious  algebraic  manipulations. 

Thus  far,  very  little  work  has  been  done  to  generalize  any  of  the 
above  results  to  any  damped  structural  member  with  arbitrary  non- 
homogeneous  boundary  conditions,  in-span  conditions,  and  non-self-adjoint 
equations.  This  is  true  even  for  viscous  damping  which  is  the  most  pop- 
ular model  because  of  its  mathematical  convenience  rather  than  its 
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physical  exactness.  O' Kelley  [7]  did  use  a general  approach  to  resolve 
the  dynamic  response  of  viscously  damped  continuous  members,  but  his 
remarks  were  restricted  to  self-adjoint  systems  with  homogeneous  boundary 
conditions  and  no  in-span  conditions.  Meirovitch  [2]  also  considered 
self-adjoint  systems  of  equations  for  both  viscous  and  structural  damping 
with  only  homogeneous  boundary  conditions  and  no  in-span  conditions.  The 
most  general  formulation  was  presented  by  Pfennigverth  [10]  and  later 
extended  by  Cinelli  and  Pilkey  [11].  However,  only  undamped  continuous 
structural  members  were  considered. 

In  this  work,  a comprehensive  theory  is  presented  for  the  dynamic 
response  of  continuous  damped  structural  members.  A general  set  of 
formulas  is  derived  that  explicitly  provides  the  ingredients  necessary 
to  form  the  modal  solution  due  to  arbitrary  loadings.  These  ingredients 
include  the  necessary  normal  modes  (shown  to  possess  the  proper  orthog- 
onality) and  the  time-dependent  uncoupled  coordinates.  It  is  shown 
that  a very  general  type  of  problem  can  be  resolved  with  this  approach. 

That  is,  linear  structural  members  can  be  self-adjoint  or  non-self-adjoint, 
possess  homogeneous  or  non-homogeneous  boundary  conditions,  or  have 
in-span  conditions  such  as  in-span  supports.  In  addition,  members 
with  two  or  three  independent  spatial  variables  may  be  treated  as  easily 
as  those  with  one  spatial  dimension.  Knowing  only  the  differential 
equations  of  the  structural  member,  the  dynamic  response  is  explicitly 
written  out,  thus  avoiding  the  tedious  algebraic  manipulations  which  are 
now  required  when  solving  for  a particular  member.  Finally,  this 
general  approach  also  exposes  the  mathematical  limitations  imposed  by 
the  physical  model  which  prohibit  the  uncoupling  of  some  damped  structural 
members. 
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General  Formulation  for  the  Dynamic  Response  of  Continuous  Structural 
Members  with  Viscous  Damping 

The  starting  point  of  this  general  formulation  is  the  governing 
differential  equation  of  motion  of  the  structure.  The  class  of  viscously 
damped  structural  members  to  be  considered  may  be  defined  by  the  following 
equation  of  motion: 


Du(x, t)  = l A. (x)3.u(x,t)-F(x,t) 
j=l  J J 

with  initial  conditions: 


(la) 


u(x,0)=Uq(x) 


3u(x,0)  = uQ(x) 


(lb) 


and  the  time-dependent  boundary  and  in-span  conditions. 


L^u(x, t) 


P.(x,t) 


on  S,  the -body  surface 
and  in-span  locations 


(lc) 


where. 


F(x,t)  = N-dimensional  column  vector  of  body  forces 
D(x),  L^(x)  = N-square  spatial  matrix  linear  differential  operators 
A. (x)  = N-square  spatial  matrix 

J c 

P^(x,t)  = N-dimensional  column  vector  of  non- homogeneous 
boundary  and  in- span  conditions 

u(x,t)  = N-dimensional  column  vector  of  dependent  state  variables 

Uq(x),Uq(x)  = N-dimensional  vectors  of  prescribed  state  variable 
initial  conditions 

3^  >=  3(j3)/3t(j) 

and  x represents  the  independent  spatial  coordinates  x3 , X2 , and  x3. 
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Since  equations  (1)  are  matrix  equations,  they  represent  N distinct 
governing  differential  equations  of  motion  and  their  associated  boundary, 
in-span,  and  initial  conditions.  Mote  that  lumping  the  boundary  and 
in-span  conditions  into  one  equation  (equation  (lc))  is  reasonable  be- 
cause boundary  terms  may  really  be  considered  in-span  conditions  at 
the  boundary  of  the  member.  Furthermore,  the  case  of  homogeneous 
boundary  and  in-span  conditions  may  be  considered  to  be  a special  case 
of  this  development  wherein  P^(x,t)  is  set  equal  to  zero. 

In  order  to  determine  the  transient  response  using  a modal  analysis, 
normal  modes  of  the  structural  member  must  be  known.  Depending  on  whether 
the  member  possesses  non-proportional,  proportional,  or  no  damping,  the 
normal  modes  used  in  the  modal  expansion  will  be  different.  It  has  been 
shown  by  Caughey  and  0' Kelley  [8]  that  damped  continuous  linear  systems 
can  be  uncoupled  by  either  undamped  or  damped  modes  and  these  constitute 
two  mutually  exclusive  classes  of  problems.  Systems  which  require  only 
undamped  modes  for  uncoupling  are  called  classical  and  all  others  are 
non-classical.  Necessary  and  sufficient  conditions  have  also  been  formu- 
lated to  cetermine  when  a viscously  damped  self-adjoint  linear  continuous 
system  is  classical  or  non-classical  [8]. 

Solution  cf  a so-called  free  vibration  problem  will  yield  the 
normal  modes  and  corresponding  eigenvalues  necessary  for  modal  analysis. 
Classical  normal  modes  will  result  if  damping  is  excluded  from  the  free 
vibration  problem;  non-classical  modes  if  damping  is  included.  3y  defi- 
nition, the  damped  or  undamped  free  vibration  problem  represents  motion 
in  which  all  external  forces  and  prescribed  deflections  have  been  set 
equal  to  zero.  When  the  prescribed  conditions  occur  on  the  boundary, 
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they  are  called  non-homogeneous  boundary  conditions.  Otherwise,  they 
constitute  in-span  conditions.  A member  undergoing  free  notion  will 
vibrate  due  to  forces  inherent  in  it,  having  been  set  in  motion  by 


prescribed  initial  conditions. 

For  an  undamped  member,  it  can  be  shown  that  the  free  motion 
response  has  a separable  form  of  solution  given  by 

iX  t 

u(x,t)  = i^n(x)e  n • (2) 

where  X is  the  natural  frequency  and  <|>  (x)  is  the  vector  of  undamped 
n n 

(or  classical)  mode  shapes.  The  undamped  governing  equation  of  motion  is 
given  by  setting  the  damping  matrix  A^  = 0 in  equation  (la) . By  substi- 
tuting the  solution  given  by  equation  (2)  into  the  governing  equation, 
the  undamped  free  vibration  problem  is  obtained 

« 

Dili  (x)  = -(X  )2  A2 (x)1^  (x)  (3a) 

n n n 

3 

with  the  boundary  and  in-span  conditions 

i 

L 4>  (x)  = 0 on  S (3b) 

l n 

It  will  be  shown  that  only  the  classical  modes  (x)  are  needed  to 

n 

express  a modal  series  solution  of  equation  (1)  for  undamped  or  pro- 
portionally damped  systems. 

Guided  by  equation  (2) , when  a member  has  non-proportional  viscous 
damping  the  damped  free  motion  solution  may  be  assumed  to  have  the  form 

iui  t -y  t 

u(x,t)  = vQ(x)e  n e n (4) 


L 
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where  e represents  the  temporal  decay  due  to  damping  in  the  struc- 
tural member.  By  defining  the  complex  frequency  to  be 


S * -Y  + iu) 
n n n 


equation  (4)  may  be  written  as 


u (x,  t ) = v (x)  e 
n 


in  analogy  with  equation  (2).  Now  and  v^Cx)  represent  the  eigenvalue 
and  vector  of  mode  shapes,  respectively,  corresponding  to  the  damped 
structural  member.  Note  that  in  contrast  to  the  undamped  member,  the 
eigenvalues  and  mode  shapes  are  comp lex- valued  for  underdamped  structural 


members. 

As  was  done  for  the  undamped  case,  the  damped  free  vibration 
problem  may  now  be  derived  by  substituting  equation  (6)  into  equation 
(la)  and  setting  all  applied  loadings  equal  to  zero.  This  results  in 
the  equation 


Dv  (x)  - l (S  )j  A.  (x) v (x) 
n . . n j n 


with  che  boundary  and  in-span  conditions 


L.v  (x)  - 0 
i n 


Equations  (3)  and  (7)  constitute  eigenvalue  problems.  The  undamped 
free  vibration  problem  led  to  equation  (3),  which  for  self-adjoint 
systems  is  a higher-order  generalization  of  the  traditional  Sturm- 
Liouville  boundary  value  problem.  When  damping  is  included  in  the  free 
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vibration  problem,  a non-classical  eigenvalue  problem  results  in  which 
the  eigenvalue  appears  non-linearly.  Regardless  of  whecher  the  system 
is  classical  or  non-classical,  it  will  be  assumed  that  the  solution  con- 
sists of  a denumerably  infinite  set  of  eigenvalues  and  eigenfunctions, 
which  are  complete.  Also  since  repeated  eigenvalues  are  of  limited 
practical  use,  only  zero  multiplicity  will  be  considered  in  this  work. 
Hence,  the  corresponding  eigenfunctions  will  be  assumed  to  be  linearly 
independent . 

The  previous  undamped  and  damped  free  vibration  problems  may  be 
combined  into  the  single  general  form: 

2 

Du  (x)  * l (-1)  ^ (u  )■*  A (x)u  (x)  (8a) 

n j-a  n j n 

with  the  boundary  and  in-span  conditions 


L.u  (x)  » 0 
l n 


on  S 


(8b) 


Here  u^  and  are  the  eigenfunction  and  eigenvalue  which  can  be  classical 
or  non-classical,  depending  on  whether  damping  is  included  in  the  free 
motion  problem.  That  is,  assign  the  following  significance: 


(u  >H_)  * 
n n 


(*?  ,X  ),  Undamped  Free  Motion 
n n 


(v  ,S  ),  Damped  Free  Motion 
n n 


(9) 


Also,  for  no  damping  set  a = 2 and  when  damping  is  present  assign  a = 1. 
The  remaining  discussion  in  this  work  will  use  this  notation  whenever 
possible  to  enhance  the  general  approach  underlying  this  work. 
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Associated  with  both  the  undamped  and  damped  eigenvalue  problems, 
an  adjoint  eigenvalue  problem  may  be  described.  Using  the  general 
notation  in  equation  (9),  the  adjoint  problem  takes  the  form 


D uq(x)  - l (-1) 
j-a 


n j n 


(10a) 


with  homogeneous  adjoint  boundary  and  in-span  conditions 


L u (x) 
i n 


where 


D(x) , Li(x) 


Aj  (x) 


u (x) 
n 


on  S (10b) 


N-square  algebraic  adjoint  spatial  matrix  linear 
differential  operators 

N-square  algebraic  adjoint  spatial  matrix 

N-dimensional  column  vector  of  algebraic  adjoint 
eigenfunctions  corresponding  to  the  nth  eigenvalue, 
either  classical  or  non-classical  (see  equation  (9)) 


The  eigenvalue  of  the  original  free  vibration  problem  is  also  the  eigen- 
value of  the  adjoint  problem,  for  both  undamped  and  damped  cases  (see 
References  (12)  and  (13)).  Note  that  in  equation  (10a)  the  matrix  A^  (x) 
has  been  restricted  not  to  be  a differential  operator;  therefore,  the 
algebraic  adjoint  (x)  is  equivalent  to  the  transpose  of  A^(x). 

The  adjoint  operator  and  the  adjoint  boundary  and  in-span  conditions 
denoted  by  D(x)  and  L^(x)  in  equation  (10a)  are  not  known  a priori  from 
the  original  statement  of  the  eigenvalue  problem.  However,  these 
operators  are  formally  shown  to  exist  by  the  use  of  various  forms  of 
the  Green's  identity.  One  form  which  is  valid  for  homogeneous  boundary 
conditions  and  no  in-span  conditions  is  given  by 
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(11) 


<u  , D u > — <u  , D u > m 0 
n m oi  n 


where  the  notation  <u  , D u > denotes  the  inner  product, 

a m 


<u  , D u > * / u • (D  u )dx 

n m , ‘ n m 

domain 


and  u and  u represent  either  the  classical  or  non-classical  vector  of 
m n 

mode  shapes  and  the  associated  vector  of  adjoint  modes,  respectively. 

When  in-span  conditions  are  present,  it  can  be  shown  [13]  that  the 
Green's  identity  takes  the  form, 

<u_»  Du>=<u,Du>+  B(u  , u ) (13) 

cl  m m n n m 

where  B(u  , u ) is  a bilinear  function  of  the  state  variable  vectors 
n m 

u^  and  uq  which  represents  the  boundary  term  and  in-span  condition  at 

each  location.  Note  that  the  eigenfunctions  u and  u must  be  differenti- 
al n 

able  to  the  extent  demanded  by  the  operators  6 and  D when  no  in-span 
conditions  are  present.  Otherwise,  these  functions  need  only  be 
differentiable  at  all  locations  not  coinciding  with  the  in-span  conditions. 

Using  the  proper  form  of  the  Green's  identity,  it  may  be  used  to 
identify  the  adjoint  free  vibration  problem  for  any  classical  or  non- 
classical  structural  member.  Briefly,  the  left-hand-side  of  equation  (13) 
is  formed,  and  then  integrated  by  parts  so  that  the  operator  6 and 
bilinear  functional  B are  identified  for  any  member.  Several  practical 
illustrations  of  deriving  equations  (10)  are  provided  in  Reference  [13]. 

Note  that  in  general,  the  solution  of  the  adjoint  problem  given  symbolically 
by  equation  (10)  will  be  different  from  that  of  the  original  free  motion 
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I 


in  equation  (8).  However,  when  the  problem  is  self-adjoint  these  two 
solutions  will  be  identical.  The  advantage  of  this  is  obvious:  only 
one  eigenvalue  problem  must  be  solved.  Moreover,  an  orthogonality 
relation  which  involves  only  the  original  eigenfunctions  is  needed  to 
uncouple  the  transient  response. 

The  modal  analysis  solution  of  a dynamic  vibration  problem  is 
characterized  by  the  form 


u(x, t)  = l q^COu^Cx)  (14) 

m*l  “ 

where  u(x,t)  is  the  vector  of  state  variables  describing  the  deflection 

and  internal  forces  present  in  the  member  at  any  given  location  and  time. 

The  normal  modes  or  eigenfunctions  u (x)  will  be  assumed  to  be  known. 

n 

In  the  classical  case,  u (x)  can  often  be  determined  analytically  in 

in 

closed  form  for  many  simple  members  [2].  However,  when  non-proportional 
damping  is  present  in  a structural  member  an  analytical  expression  for 
the  eigenfunctions  cannot  be  found  and  instead  a numerical  technique  must 
be  used.  When  the  structural  member  has  one  spatial  direction,  transfer 
matrices  are  particularly  well-suited  to  providing  an  accurate  means  of 
solution.  A description  of  this  technique  for  numerically  evaluating 
the  eigenfunctions  will  be  given  in  a later  section. 

In  this  work,  it  is  assumed  that  the  modal  series  solution  given  by 
equation  (14)  is  convergent  for  all  physically  admissible  problems. 
Furthermore,  it  will  be  assumed  that  the  series  is  uniformly  convergent 
except  at  locations  where  discontinuities  exist  such  as  in-span  supports. 
At  these  locations,  the  series  may  or  may  not  converge.  Hence,  it  may 
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be  concluded  that  the  dynamic  response  will  be  determined  if  the  proper 
eigenfunctions  are  shown  to  be  theoretically  orthogonal,  and  the  so-called 
normal  coordinates  q^(t)  can  be  found  via  an  uncoupled  differential 
equation  for  each  m-coordinate. 

To  find  the  dynamic  solution  of  a general  damped  structural  member 

using  modal  analysis,  an  orthogonality  relation  must  be  shown  to  exist 

between  the  eigenfunctions.  Because  both  the  classical  and  non-classical 

free  vibration  problems  can  be  non-self-adjoint , orthogonality  will  be 

shown  to  prevail  between  the  vector  of  original  eigenfunctions  u^tx) 

and  the  associated  vector  of  adjoint  functions  u (x) . This  is 

in 

commonly  referred  to  as  a biorthogonality  relation.  When  the  free  motion 

problem  is  self-adjoint,  the  biorthogonal  form  reduces  to  an  orthogonality 

condition  between  the  original  eigenfunctions  only.  In  this  work  both 

forms  will  be  referred  to  as  an  orthogonality  relation. 

To  derive  the  general  orthogonality  relation,  begin  with  the  form 

of  Green's  identity  given  in  equation  (11).  Note  that  even  if  the  member 

possesses  in-span  conditions,  the  extended  Green's  identity  in  equation 

% 

(13)  reduces  to  the  homogeneous  form  by  Btu^.u^)  = 0.  Substitute  equations 
(8a)  and  (10a)  into  equation  (11)  to  give 


3-  » 


(-1)  ^a~1^  C(u_)^  < u,Au>  - (u  )-*  < u , A . u > ] - 0 
njm  n mjn 


m 


(15) 


Another  form  of  the  Green's  identity  may  be  given  by 


<u  ,A,u  > <u  ,A.u  > j ■ 1,2 

mjn  n j m J ’ 


(16) 
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which  is  valid  for  A^ , a matrix  of  scalar  spatial  elements.  A common 
inner  product  may  now  be  factored  out  of  equation  (15)  by  applying 
equation  (16),  so  that  equation  (15)  becomes 


I 

j-a 


(-1)(C“1)[(  u)J  - (U  )J]  <u,A  u > 
m n n j m 


(17) 


The  form  of  the  orthogonal  relation  depends  on  whether  the  system  is 
classical  or  non-classical.  Considering  the  former  case,  let  a = 2 and 
u^(x)  and  take  on  the  values  given  by  equation  (9).  Then  equation 
(17)  may  be  rewritten  as 


(g  -6  ) <<p  ,A2i P > =0 
n m n m 


(13) 


where  3 = X2.  When  3 t4  8 the  inner  product  must  be  zero  to  satisfv 

n n n m 

equation  (18);  otherwise,  it  does  not  vanish.  Thus,  the  orthogonal 
relation  for  classical  systems  may  be  expressed  in  the  form: 


«l/  ,A2<I>  > 3 6 Q 

n m mn  n 


(19) 


where  represents  the  classical  norm  given  by 


(20) 


and  5 is  the  Kronecker  delta.  Note  that  in  the  unlikely  event  that 
mn 

the  norm  is  zero,  the  vector  of  adjoint  eigenfunctions  may  be  replaced 

by  its  complex  conjugate.  The  functions  and  ’Jj  would  then  be  orthogonal 

n m 

in  the  so-called  Hermitian  sense. 

For  non-classical  systems,  a = 1 and  u and  u become  v and 

n n n 

Sn,  the  vector  of  damped  eigenfunctions  and  correspnding  complex 
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eigenvalues  in  equation  (15).  After  applying  the  same  argument  that  led 
to  equation  (19),  non-classical  orthogonality  may  be  expressed  by 


(S  -S  ) [(S  +S  ) <v  ,A?v  > + <v  , AiV 
mn  mn  n 4 m n 1 m 


5 N 
mn  n 


(21) 


where  now  N 

n 


is  the  non-classical  norm  given  by 


N * 2S  <v  ,A?v  > + <v  ,Aiv  > 
n nn  n n n 


(22) 


Mote  that  Aj  and  may  contain  complex  elements  and  therefore,  both 
norms,  and  Q^,  may  be  complex- valued. 

To  complete  the  modal  solution  as  given  by  equation  (14)  the 
uncoupled  coordinates  9m(b)  must  be  determined.  This  is  accomplished  by 
transforming  the  governing  equation  of  motion  into  an  uncoupled  differ- 
ential equation.  Uncoupled  here  means  that  a different  equation  exists 
for  each  m-eigenvalue,  independent  of  all  ether  values.  To  effect  this 
transformation,  begin  with  the  Green's  identity  in  the  general  extended 
form 


<u  ,Du>  - <u,D  u > * 3 (u.  ,u)  (23) 

m mm 

This  equation  is  valid  for  a member  with  non-homogeneous  boundary  conditions 
and  in-span  conditions  (B  f 0)  as  well  as  for  homogeneous  boundary  and  in- 
span conditions  (B  = 0).  Substitute  equation  (la)  and  (10a)  into  the  above 
equation  to  obtain 


2 ' 2 

l 3,  <u  ,A  u>  " l (-D  (a_1)  (u) ^ <u,A.u  > - 
|-1  j m 3 j-a  m 3 m 


<u  ,F>  + B ( u ,u) 
m m 


(24) 
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The  bilinear  form,  once  identified  for  a particular  example,  will  retain 
its  form.  By  applying  the  following  form  of  the  Green's  identity, 


<u,A.u  > * <u  ,A  u> 
j m m j 


j **  1,2 


equation  (24)  may  be  rewritten  as 


l 3.  <uA  u>  - ,1  (-l)(a-1) (y  <u  ,A  u> 
|=1  i m 2 a m m 2 


<u  ,F>  + B(u  ,u) 
m m 


To  reduce  this  equation  further,  the  value  of  a must  be  known.  For  a 

classical  system,  a = 2 and  u and  y are  replaced  by  the  adjoint  classi- 

m m 

cal  mode  shape  <(>  and  the  corresponding  eigenvalue  X . Noting  these 
m m 

changes  and  expanding  the  sum,  equation  (26)  becomes 


3 2 + (X  ) j <;fc  ,Aou>  + 3<’^  ,Aiu> 
mm  m 


<i0  ,F>  + 30l>  ,u) 
m m 


A common  inner  product  may  not  be  factored  out  due  to  the  presence  of 
the  damping  term  A^ . At  this  point  the  proportionality  assumption  must 
be  invoked. 


A^u  * aA£U  - b D u 


where  a and  b are  constants  of  proportionality.  That  is,  the  damping 
terms  (A^)  are  assumed  to  be  proportional  to  the  mass  (A2)  and/or 


F 


stiffness  (D)  of  the  member.  Another  form  of  this  condition  may  be 
found  by  multiplying  equation  (28)  by  and  integrating  over  the 
domain  to  obtain 


<*m,Aiu>  - a <tf»m>A2u>  - b <^n,Du>  (2 

Apply  the  extended  Green's  identity  to  the  last  term  on  the  right-hand 
side  and  then  substitute  the  classical  form  of  equation  (10a)  into 
equation  (29).  After  some  rearrangement  obtain  the  equation 


<i  ,Aju>  = (a+bX2)  <’p  ,A2u> 
m m m * 

- b B(5  ,u) 
m 

With  this  form  of  the  proportionality  assumption,  return  to  equation 
(27).  After  substituting  equation  (30),  obtain  the  equation 


£ (t)  + (a+bX2) 
in  m 


m 


+ X 


2r 

m’m 


“ + (l+b3)B(^m,u) 

- F 

m 

where  ( ) * 3 ( )/3t  and  where  has  been  assigned  the  value 

Cm(t)  " <i^b»A2u>  (32) 

Note  that  is  a function  of  time  alone  because  the  spatial  variables 
have  been  eliminated  by  the  definite  integral.  The  solution  of  equation 
(31)  may  be  easily  found  to  be 
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-A  C t Ac 

c (t)  * e m m [cos  a t + -2-^  sin  a t]  5 (0) 
m m mm 

m 


-Act  sin  a t 

+ e m m [ ] I (0) 

ci  in 

m 


(t) 


m is 


sin  a (t-t) 
m 


a 

in 


dt 


(33a) 


where 


C 

m 


- (~ 
2 ^A 


m 


+ b A ) 
a 


(33b) 


and 


a 

m 


A 

a 


(33c) 


where  the  initial  conditions,  5 (0)  are  given  by  equation  (32)  by  re- 
placing  u by  and  u^,  respectively.  Note  that  as  a special  case  of  a 
classical  system,  equation  (31)  may  be  used  to  represent  an  undamped 
member.  For  this  case,  set  a = b = 0 in  equations  (31)  and  (33)  and 
continue  to  use  equation  (32). 

In  an  identical  manner,  an  expression  for  the  normal  coordinates 
for  a non-classical  system  can  be  derived.  Beginning  with  the  general 
equation  (26),  let  a * 1 and  u^  and  y^  then  become  the  adjoint  non- 
classical  mode  shape  and  the  complex  eigenvalue,  respectively.  When 
these  changes  are  introduced  into  equation  (26)  and  the  sums  are 
expanded,  it  may  be  expressed  as 
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(34) 


(3-S  ) [ <v  ,Aiu>  + (3+S  ) <v  ,A2u>  ] 
m m mm 

= <v  ,F>  + B(v  ,u) 
m m 


The  term  in  square  brackets  is  a function  of  time  alone  so  that  it  may 
be  represented  by  the  variable 

D (t)  = <v  ,Aiu>  + <v  ,A23v>  + S <v  ,A2v>  (35) 

m *tii  m mm 

Using  the  above  definition,  equation  (34)  may  be  written  as 


n (t)  - S n = <v  ,F>  + B(v  ,u)  (36) 

m m in  m m 

This  equation  is  essentially  an  ordinary,  first  order,  uncoupled  dif- 
ferential equation  which  has  the  solution 


n (t) 
m 


S t 
m 
e 


n (0)  + 

m 


Sm(t-T)  [ <v  ,F  + B (v  >u)]  dr  (37) 

J e m m 

0 


and  n (0)  is  found  bv  evaluating  equation  (35)  at  t = 0. 
n 

The  solutions  given  by  equations  (33)  and  (37)  represent  the 
first  step  in  finding  the  uncoupling  coordinates  ^(t)  necessary  to  com- 
plete the  modal  solution  given  by  equation  (14).  To  avoid  confusion, 

let  q (t)  assume  the  values 
m 


q_(t) 

m 


fm(t> 

®m(t) 


Classical  Systems 
Non-Classical  Systems 


(38) 


18 


I 

I 

I 

I 

I 

I 


I 

! 




i »■«■»  i.W»nrjWW3lBHjl 


in  equation  (14).  Considering  first  the  non-classical  systems,  recall 
that  the  ultimate  goal  is  to  find  g (t)  in  the  series  solution, 

HI 

00 

u(x,t)  « l gm(t)vm(x)  (39) 

m*l 

Since  both  n and  g are  functions  of  time  alone,  a relationship 
m m 

between  them  is  sought.  For  no  in-span  conditions  equation  (39)  is 
assumed  to  be  uniformly  convergent.  When  in-span  conditions  are  present, 
discontinuities  of  some  of  the  state  variables  will  occur  at  these  loca- 
tions. However,  at  all  other  points  along  the  member,  the  modal  expan- 
sion may  be  assumed  to  remain  uniformly  convergent.  That  is,  the  state 
variables  will  be  assumed  to  be  continuous  throughout  all  subintervals 
which  are  punctuated  by  these  discontinuities.  Such  functions  are  often 
referred  to  as  sectionally  continuous.  For  any  well-defined  physical 
problem,  this  modal  expansion  may  converge  very  slowly  or  even  diverge 
at  these  points  of  discontinuity.  However,  the  analyst  is  frequently 
not  interested  in  the  response  at  precisely  these  points  so  that  the 
formulation  in  equation  (39)  is  still  useful  as  long  as  these  loca- 
tions are  excluded  from  the  analysis.  Then,  at  all  points  for  which 
the  series  is  uniformly  convergent,  it  may  be  differentiated  with  respect 

to  time.  Guided  by  equation  (37),  g (t)  will  have  an  exponential  form 

m 

so  it  follows  that 


” l smgm(t)vn,(x) 

ot  , m m m 

m=l 


(40) 


Substitute  equations  (39)  -ind  (40)  into  the  definition  of  n (t) , 
then  by  assuming  that  the  integral  of  an  infinite  series  is  equivalent 
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to  an  infinite  series  of  the  integrals,  this  equation  may  be  rewritten 
as 


n (t)  - 7 g [<v  , AiV  > + S <v  ,A 2v  > 
nnv  ' ^ n 1 m m n m 

m~l 


+ S <v  ,A2v  >] 
n n m 


(41) 


Using  non-classical  orthogonality  and  the  associated  norm  in  equations 
(21)  and  (22) , equation  (41)  becomes 


g_(t) 

in 


n (t) 
m 

N 

m 


(42) 


Since  n (t)  is  known  from  equation  (37),  g (t)  is  also  known.  Hence, 
m m 

the  dynamic  response  may  be  written  out  as  shown  in  equation  (39). 

For  the  classical  system,  the  same  procedure  may  be  followed  to 


give 


fm(t) 

in 


*m(t) 

m 


(43) 


a.  well-known  result  for  undamped  structural  members.  The  modal  solution 
then  takes  the  form 


00 

u(x,t)  = l ~Q * m(x) 

m=l  ^m 

In  summary,  the  dynamic  response  of  any  structural  member  with 
equations  of  motion  given  by  equation  (1)  may  be  found  by  applying 
either  of  the  following  two  series: 
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’ I 


u(x,t) 


" ^ (t) 

7 — r (x)  , Classical 

•,  y m 

m=l  m 


00  n (t) 

Z -i — v (x) 
, N m 
m*l  m 


Non-Classical 


(45a) 


(45b) 


where  C (t)  and  n (t)  are  given  by  equations  (33)  and  (37),  respectively, 
mm 


Acceleration  Method  of  the  Dynamic  Response  of  Continuous  Structural 
Members  with  Viscous  Damping 

In  the  previous  section,  the  dynamic  response  of  a completely 
general  viscously  damped  structural  member  was  formulated  as  an 
infinite  series  solution.  This  form  of  the  solution  is  often  referred 
to  as  the  displacement  method.  The  advantage  of  this  approach  lies  in 
that  it  will  usually  converge  to  the  response  with  only  a small  number 
of  terms.  If  the  series  solution  requires  many  terms,  it  is  often  pos- 
sible to  overcome  the  slow  convergence  by  modifying  the  displacement 
solution.  One  technique  which  may  be  used  to  achieve  this  goal  is 
called  the  acceleration  method.  In  this  section,  this  form  of  the 
modal  solution  will  be  derived  for  a structural  member  with  proportional, 
non-proportional,  or  no  viscous  damping. 


Acceleration  Method  for  Proportionally  Damped  Structural  Members 

The  structural  member  treated  in  this  subsection  possesses  the 
so-called  proportional  damping,  which  satisfies  the  condition  given  in 
equation  (30).  When  a * b = 0 the  member  is  undamped,  and  so  undamped 
members  may  be  considered  to  be  a special  case  of  proportional  damping. 
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As  mentioned  in  the  previous  section,  proportionally  damped 
structural  members  are  categorized  as  classical  systems.  Hence,  only 
the  undamped  modes  are  required  to  uncouple  the  dynamic  response  which 
is  given  by  equation  (44).  Recall  that  the  uncoupling  coordinates  £m(t) 
to  obtain 

c (t)  - - — [5  + (a+bX  ) | - Fj  (46) 

m ^2  m m m m 

m 


Substitute  this  equation  into  equation  (44)  and  arrive  at  the  expression 


4. 

® tp  F‘  ® i(i 

u(x,t)  * l -2-2-  - l — 2_  [5  + (a+X2b)  \ ] 


m=l  X20  m=l  X2Q 
in  mm 


(47) 


The  first  term  on  the  right-hand  side  of  this  -equation  is  a function 
of  the  spatial  and  temporal  variables.  However,  time  only  enters 

4. 

through  the  generalized  forcing  function  F'  which  represents  applied 

m 

loadings  and  non-homogeneous  boundary  and  in-span  conditions.  Because 

4a 

these  are  prescribed,  F^  may  be  determined  at  every  instant  of  time 
for  which  this  first  term  may  be  considered  to  be  a function'  of  x 
alone.  With  this  understanding,  define  this  expression  as  a static  or 
quasi-static  term,  given  by 


l 

m-1 


i|»  F* 
m m 

X2Q 

m m 


(48) 


For  additional  convenience,  define 


IQ 


+ (a+bX2)  k ] 
m m 


(49) 
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I 

I 

I 

I 

I 


so  that  Equation  (47)  may  be  rewritten  as 


C+(t) 


i(x,t)  - v‘  + l (x) 

s i Q_  ni 
m-1  th 


Note  that  u(x,t),  v , and  (x)  are  column  vectors  representing  the 

s m 

dynamic  static,  and  classical  mode  state  variables  of  the  member, 
respectively. 

The  underlying  goal  of  the  acceleration  method  is  to  extract  a 
static  solution  from  the  modal  expansion.  In  effect,  the  acceleration 
solution  achieves  a "jump"  on  the  dynamic  response  by  beginning  with 
the  static  solution.  To  define  the  terms  in  equation  (50)  in  greater 
detail,  a specific  solution  form  is  needed.  Recall  the  solution  for 

E (t)  given  by  equation  (33).  To  extract  the  static  solution,  integrate 

m 

the  last  term  in  equation  (33a)  by  parts  twice.  Substitute  the  resulting 
equation  into  equation  (45a),  and  then  compare  this  expression  with 
equation  (50)  for  determining  u(x,t).  In  this  way,  the  alternative 

M 

expression  for  Cm(t)  may  be  identified  as 


^(t)  ct  c*  r 

C1  (t)  ■ — + e [cos  a t + sin  a t]  [£  (0)  - — 

m (1-?2)X2  m am  mm  a2 


F ‘ (0) 


-X  ? t sin  a t , dF  (0) 

+ e " “ — -S-  [I  CO)  - i-  -B—  I + 

am  7 at 

m a.*- 

m 


- ~ f ([-7+  U + b\2)^-+  ‘(\  r )2]  F'(T) 
2 Jj2  max  mm  m 


a2  i dx2 
m 0 


mm  m 


-X  ? (t-x)  sin  a (t-x) 
mm  m . 

e — ; dT 

m 
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To  summarize,  an  alternative  solution  of  the  transient  response  problem 

has  been  derived  in  equation  (50).  The  static  solution,  v"*”  is  obtained 

s 

from  equation  (48)  and  the  temporal  coefficient  C+(t)  is  given  in 

m 

equation  (51).  Note  that  v’  may  also  be  found  by  solving  equations 

s 

(1)  in  which  u =*  v^  and  k\  * A2  = 0 by  any  suitable  method. 

As  special  cases  of  the  above  acceleration  solution,  note  that  if 

the  member  was  undamped  merely  setting  ? equal  to  zero  (a  = b ■ 0)  and 

m 

a * X would  allow  uncoupling  of  the  response.  In  addition,  the  case 
m m 

of  the  homogeneous  boundary  conditions  and  no  in-span  conditions  would 
affect  only  the  value  of  the  generalized  force  term  defined  in  equation 
(31).  That  is,  set  B(<Ji  ,u)  equal  to  zero. 


Acceleration  Method  for  Non-Proportionally  Damped  Structural  Members 
The  assumption  of  proportional  damping  is  relaxed  in  this  sub- 
section. Hence,  the  dynamic  response  of  such  structural  members  will 
be  uncoupled  with  the  non-classical  or  damped  mode  shapes.  Recall 

that  for  general  viscous  damping  the  uncoupling  coordinates  n (t)  are 

in 

found  by  solving  equation  (36).  To  find  the  acceleration  method  for 
this  case,  begin  by  solving  this  equation  for  n (t): 


n (t) 

m 


F ) 

m 


(52) 


Where  the  non-classical  form  of  the  generalized  forcing  function, 


F 

m 


(t) 


< ^ ,F  > + »*0 

in  m 


(53) 


has  been  employed.  Upon  substitution  of  equation  (52)  into  equation 
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(45b) , obtain  the  expression  for  the  dynamic  response 

• 

00  V F 00  n 

/ \ r mm.  p m 

u(x,t)  ■ -I  s N l s N vm 

m=l  m m m*l  m m 


(54) 


In  identical  fashion  to  the  previous  subsection,  the  first  series  in 
this  equation  may  be  recognized  as  being  quasi-static  so  that  it  may 
be  rewritten  in  the  form 


u(x,t)  “ v + l 


C (t) 


m 


V v 

, N m 
m*l  m 


(55) 


where  v is  the  static  solution  for  the  non-classical  system, 
s 

Following  the  pattern  of  the  previous  subsection,  the  solution 
for  the  temporal  coefficient  n (t)  in  equation  (37)  is  integrated  bv 

til 

parts  twice  to  extract  the  static  portion  of  the  response.  Then  the 

resulting  equation  for  n (t)  is  substituted  back  into  the  displacement 

m 

modal  expansion  in  equation  (45b).  This  equation  is  then  compared  with 
equation  (55),  the  static  portion  of  the  solution  is  identified,  and 
the  coefficient  C (t)  is  given  by 

TO 


. dF  (t)  S_t 
C (t)  ■ - — ” + e 


dF  (0) 


m 


;2  dt 


m 


mC  fl  dV'- 

Ls*  dt 

m 


Fm(0)  "I 

I— + v°  n 

m J 


(56) 


ft  Sm(t_T)  d2F 
e m 


dr 


m 


dx^ 
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Hence,  the  acceleration  form  of  the  modal  solution  is  given  by  equation 


(55)  where  C (t)  is  found  from  the  above  eauation.  No  simplifying 
m 

assumptions  have  been  made  in  this  derivation,  so  that  the  structural 
member  can  possess  general  (non-proportional)  viscous  properties,  have 
time-dependent  non-homo geneous  boundary  and  in-span  conditions,  and 
accept  any  arbitrary  loadings.  For  a structural  member  f^r  which  the 
usual  modal  expansion  is  prohibitive  because  of  the  many  terms  needed 
for  adequate  convergence,  this  acceleration  form  of  the  solution  may 
still  provide  a viable  means  of  finding  the  dynamic  response  using 
modal  analysis. 

Dynamic  Response  of  a Viscoelastic  Structural  Member  on  a Proportional 
Viscous  Foundation:  Part  I 

In  this  section  the  response  of  a general  viscoelastic  structural 
member  will  be  found  using  modal  analysis.  The  constitutive  relation 
to  be  used  will  be  of  the  hereditary  integral  variety.  However,  it  can 
be  shown  that  the  final  solution  may  be  altered  to  allow  a differential 
operator  form  of  the  material  law  to  also  be  used  (see  Reference  [13]). 
As  mentioned  previously,  Valanis  [3]  was  the  first  to  resolve  the 
solution  of  such  a general  "viscoelasto-kinetic"  problem  into  a 
superposition  of  a viscoelastic  quasi-static  and  a dynamic  elastic 
solution.  This  section  will  go  beyond  the  work  of  Valanis  and  others 
(Robertson,  [4], [5]  etc.)  by  considering  a more  general  approach 
applicable  to  both  self-adjoint  and  non-self-adjoint  systems,  which 
may  also  possess  a proportionally  damped  viscous  foundation  (hereafter 
a viscoelastic  member  on  a viscous  foundation  will  be  referred  to  as  a 

L 
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in  this  section  that  this  restriction  is  necessary  to  uncouple  the 
dynamic  response  of  viscoelastic  structural  members  using  the  general 
approach  involving  the  Green's  identity. 

The  starting  point  of  this  general  formulation  will  again  be  the 
differential  equations  governing  the  motion  of  the  member  given  by 

2 

G*Du(x,t)  = T A (x)  3.u(x,t)  - F(x,t)  (57a) 

j-1  * J 

with  the  initial  conditions, 

u(x,0)  - Uq(x) ; 3u(x, 0)  = uQ(x)  ‘ (57b) 

and  the  time-dependent  boundary  and  in-span  conditions, 

G*L^u(x, t)  - Pi(x,t)  on  S (57c) 

where  P^(x,t)  is  a vector  containing  the  prescribed  non-homogeneous 
surface  tractions  and  homogeneous  prescribed  displacements.  The 
notation  used  for  viscous  damping  remains  unchanged  for  the  general 
viscoelastic  problem.  Hence,  proportional  viscous  damping,  such  as 
found  in  a foundation,  is  included  in  the  matrix  (x) . In  addition, 


In  fact,  it  can  be  shown  [13]  if  the  viscous  foundation  is  non- 
proportional, a modal  solution  is  not  possible. 
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equation  (57c)  contains  all  boundary  conditions  including  the  non- 
homogeneous  displacements.  Although  the  prescribed  displacements  do 
not  involve  G*,  equation  (57c)  can  be  used  to  represent  these  boundary 
conditions  because  they  are  homogeneous  and  G is  a monotonicaliy 
decreasing  function  of  time.  Note  that  equations  (57)  apply  only  to 
a linear  isotropic  viscoelastic  material,  deforming  under  isothermal 
conditions  with  a constant  Poisson's  Ratio,  and  which  is  incompressible. 
Nevertheless,  this  general  form  of  the  equations  of  motion  describe  a 
wide  class  of  viscoelastic  members  on  a proportionally  damped  viscous 
foundation.  An  example  of  such  a structural  member  will  be  shown  to 
fit  this  general  theory  in  the  Applications  section. 

Recall  that  the  solution  of  a general  self-adjoint  viscoelastic 
member  has  been  given  by  Valanis  to  be  a superposition  of  a quasi- 
static viscoelastic  and  a dynamic  elastic  solution.  Guided  by  this 
observation  and  the  fact  that  the  viscous  damping  is  restricted  to  be 
proportional,  it  is  presumed  that  the  undamped  mode  shapes  may  be  used 
to  uncouple  the  response  of  any  complicated  member.  Recall  from  a 
previous  section  that  the  classical  eigenfunctions  were  found  by  con- 
sidering the  undamped  free  vibration  problem  given  by  equations  (3).  In 
addition,  the  viscoelastic  member  may  be  non-self-adjoint,  so  that  a 
biorthogonality  relation  will  be  needed  in  deriving  the  elastic  modal 
expansion.  Hence,  use  will  be  made  of  the  adjoint  undamped  free  motion 
problem  given  by  equations  (10)  along  with  equation  (9).  Because  the 
elastic  undamped  free  vibration  problems  are  identical  to  the  proportional 
viscous  damping  case,  the  orthogonality  relation  given  by  equations 
(19)  and  (20  is  also  unchanged. 


Assume  that  the  solution  of  the  general  viscoelastic  dynamic 


I 


I 


I 

! 


problem  may  be  expressed  in  the  usual  modal  series  form. 


u(x,t)  - l h (t)  'J'  (x)  (58) 

, m m 
m*l 

The  classical  eigenfunctions,  (x) , can  be  easily  found  for  a variety 

E 

of  structural  members.  Thus,  only  the  temporal  coefficients  h (t)  need 

in 

to  be  determined  to  complete  the  modal  solution.  In  order  to  derive 
an  uncoupled  equation  for  these  coefficients,  proceed  as  follows.  The 
Green's  identity,  in  extended  form,  may  be  written  as 


< rj)  ,D(G*u)  > - < (G*u) , Dip  > =•  B(\p  , (G*u) ) 
m mm 


(59) 


The  term,  (G*u) , is  really  a function  of  both  the  spatial  and  temporal 

variables  and  its  use  in  equation  (59)  is  valid  if  it  is  differentiable 

to  the  extent  demanded  by  the  operator  D.  This  is  assumed  to  be  the 

case.  The  bilinear  form  B(^  ,(G*u))  is  determined  in  exactly  the  same 

m 

way  as  for  the  elastic  member.  That  is,  the  inner  product  < <l>  ,D(G*u)  > 

m 

is  formed  and  then  integrated  by  parts  with  respect  to  the  spatial 
coordinates.  Note  that  the  appearance  of  G does  not  effect  the  inte- 
gration by  parts  because  G is  a function  of  time.  The  boundary  and 
in-span  conditions  are  then  grouped  to  form  the  bilinear  functional 
B(^m, (G*u)) . Since  the  prescribed  displacements  on  the  boundary  are 
zero,  this  bilinear  form  contains  only  the  non-homogeneous  surfacq 
tractions  which  involve  the  relaxation  modulus  G(t).  In  fact,  it  is 
this  form  of  B(^m> (G*u))  that  prohibits  non-homogeneous  displacements 
from  being  prescribed  on  the  boundaries  because  these  may  not  be  expressed 
in  the  form  G*u. 
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Before  returning  to  equation  (59),  note  that  a property  of  the 
convolution  is  that 

G*(Du(x,t))  » (Du(x,t))*G  (60) 

where  Du(x,t)  may  be  treated  as  an  arbitrary  function  of  x and  t. 

Since  D is  a spatial  differential  operator  it  does  not  effect  G,  a 
function  of  time,  so  that 

G*Du  = Du*G  = D(u*G)  = D(G*u)  (61) 

With  this  equation,  the  governing  equation  of  motion,  equation  (57a) 
becomes 


D(G*u  = l A (x)  3 u(x,t)  - F(x,t) 

j-1  J J 

Now  substitute  equations  (10a)  and  (62)  into  the  Green’s 
equation  (59)  to  obtain 

3 < 4*  ,Ai  u > + 3o  < >Aou  > + < (G*u) , A?P  > 

m 1 * m * m ^ m 

- < -J)  ,F  > + B(t^  , (G*u) ) 

m m 

Noting  the  homogeneous  form  of  the  Green's  identity, 

< (G*u) , A24,m  > * < i(im,A2(G*u)  > 

and  that  G is  a function  of  time,  it  follows  that 


(62) 

identity  in 


(63) 


(64) 


< (G*u)  ,A2^  > » G*  < \p  yAyu  > 

m m 


(65) 
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1 

I 

Using  equation  (65)  in  equation  (63)  gives 

• 3^  < '!>  , A?u  > + 3 < i(i  ,AiU  > + (X  )^  G*  < <p  jAtU  > 

in  m m m 

■ < i|»  ,F  > + B(t{i  ,(G*u))  (66) 

m m 


by  using  the  proportionality  condition  in  equation  (30) , a common 
inner  product  may  be  factored  out  of  equation  (66)  to  give 


£ (t)  + (a  + bA*)  | + (X  )2  G*E 

m m m in  th 

= < ,F  > + B(-J>  , (G*u) ) + b8B(’i»  ,u) 
mm  m 


(67) 


where  £ (t)  has  been  defined  as  before  in  the  viscous  damping  case  in 
m 

equation  (32).  Likewise,  the  generalized  initial  conditions,  £^(0)  and 

5 (0)  are  found  from  this  equation  at  t = 0.  The  coefficients  h (t) 
m m 

are  found  as  before  for  undamped  members  to  be 


h (t) 
m 


^m(t) 

Qm 


(68) 


where  5 (t)  is  found  from  equation  (67).  Note  that  although  5 (0)  has 
m m 

the  same  form  for  elastic  and  viscoelastic  members,  the  solution  for 

£ (t)  will  in  general  be  different,  due  to  the  G*5  term  in  equation 
m m 

(67).  These  solutions  coincide  when  G*  reduces  to  unity,  i.e.  the 
viscoelastic  member  becomes  elastic. 

Although  a solution  of  the  viscoelastic  dynamic  problem  may  in 
principle  be  given  by  equation  (58) , this  form  of  the  modal  solution 
is  seldom  used  because  of  its  poor  convergence  qualities.  In  fact,  most 
authors  do  not  even  consider  the  displacement  form  as  a viable  approach. 
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Therefore,  attention  will  be  focused  on  the  more  rapidly  converging 

acceleration  (or  Williams)  method  for  a general  viscoelastic  structural 
member. 

The  acceleration  form  of  the  modal  solution  may  be  derived  by 
following  the  steps  used  in  the  previous  section.  Begin  by  rearranging 
equation  (67)  into  the  form, 


G*5 


m 


(X  ) 

m 


7 [ 5 + (a+bA2)  } - H ] 

2 m m m m 


where  is  the  generalized  force  term  given  by 


(69) 


H (t)  = 
m 


^mF  dx  + B0f>m(G*u))  + b 3(*m,u) 


(70) 


Taking  the  convolution  of  both  sides  of  equation  (53)  and  using  equation 
(68)  gives 


G*u  - £ 


- (G*5  ) 


m-1  Qm 


<1*  (x) 


(71) 


assuming  the  summation  and  convolution  may  be  reversed.  Now  introduce 
equation  (69)  into  (71)  and  after  some  rearrangement  obtain 

4>  H 


G*u  - l 

nr*l  Am  m 


* 


Z p 

m-1  A2Q  L 
mm 


5 + (a+bA2)  | 

m m m 


(72) 
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■ft  ummi  jin 


The  first  tern  in  equation  (72)  is  the  static  contribution  to  the  vis- 
coelastic solution.  As  Valanis  [3]  and  Robertson  and  Thomas  [4]  have 
pointed  out,  a dynamic  viscoelastic  problem  has  a solution  of  the  form 


00  E (t) 

u(x,t)  = vg(x,t)  + l -q ^(x)  ( 

m=l  m 

where  vg  is  the  quasi-static  viscoelastic  solution.  Equations  (72) 
and  (73)  are  equivalent  by  letting 


« 'll  H 

G*v  = £ -HUE 

3 m-1  A2Q 

m m 


oo  r. 

= - y — 5 + (a+bA2)  Z 

, , 2 m mm 

m=l 

m *—  • — 


Then  equation  (72)  becomes 


G*u  - G*v  + 7 (G*E  ) 

c “ m 0 


which  implies  equation  (73).  A more  thorough  understanding  of  the  terms 
in  equation  (76)  may  now  be  undertaken.  The  Green's  identity  may  be 
expressed  in  the  form 


< ip  , D(G*v  ) > » < (G*v  ),  D<£  > + B(-£  , (G*v  ))  (77) 

ns  s n n s 


Substitute  equations  (10a)  and  (74)  into  this  equation  to  obtain 


< <p  , D(G*v  ) 
n s 


« r-  _ __  - 

> = - y —2 — H (t)  < tp  Aoip> 

m-1  X2Q  m m’  n 

m m l—  J 


+ B(ip  ,(G*v  )) 
n s 


(78) 


where  integration  of  the  series  is  assumed  to  be  equal  to  a series  of 
the  integrals.  Using  the  orthogonality  relation  for  undamped  systems 
and  in  view  of  equations  (61)  and  (70),  equation  (78)  becomes 


< tp  , G*Dv  > = - < ip  , F > - R('|>  , (G*u) ) 
ns  n n 

- b 3 BOP  , u)  - B*0P  , (G*v  )) 
n ns 


(79) 


The  bilinear  forms  in  the  above  equation  contain  contributions  to  the 

generalized  forcing  function  from  the  non-homogeneous  tractions  on  the 

surface.  In  particular,  the  form  B(ip  , (G*u))  contains  the  non- 

n 

homogeneous  surface  tractions  as  given  by  equation  (57c-).  Therefore, 
B(tp^,  (G*vg))  implies  that  the  same  boundary  and  in-span  conditions 
apply  to  vs  when  damping  is  absent  (b  = 0) . In  addition,  a differential 
equation  valid  throughout  the  volume  is  implied  by  the  remaining  two 
terms  in  equation  (79).  That  is,  the  following  boundary  value  problem 
may  be  reduced  from  equation  (79) : 

G*D(vg (x, t) ) - - F(x,t)  (80) 

with  boundary  and  in-span  conditions 

G*LiVs(x,t)  - Pi(x,t)  on  S (81) 
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which  when  solved  will  yield  the  quasi-static  response  of  a viscoelastic 
member  used  to  obtain  the  complete  dynamic  response  as  shown  in  equation 
(73).  Note  that  this  solution  is  quasi-static  in  the  sense  that  vg  is 
a function  of  x alone  which  is  to  be  found  for  each  instant  of  time. 

To  complete  the  acceleration  form  of  the  solution,  an  uncoupled 


equation  for  E (t)  must  be  found.  Substitute  the  series  representation 
m 

in  equation  (73)  into  the  governing  equation  of  motion,  equation  (57a). 
This  gives 


G*Dv  + G*D  l E ip 
s m m 

m*l 


+ A 7 E + A v + Ao  -F 

“ tt»  m o L*  m m 


where  it  is  assumed  that  the  infinite  series  converges  at  all  points  of 
interest.  The  first  term  in  equation  (82)  is  equal  to  -F  which  can 


be  cancelled  from  either  side  of  this  equation.  Premultiply  by 
then  integrate  over  the  volume  to  obtain 


a,vc  > + i £ < * , ax*  > + < , a2; 


+ l E_  < 4>_,  A2<(i  > * l (G*E  ) < ip  , D«  > 

m-1  m n 11  m*l  m n ” 


Note  that  the  validity  of  reversing  the  convolution  and  the  summation 
sign  has  been  assumed.  Apply  equation  (3)  to  the  above  equation  and 


then  use  the  orthogonality  relation  for  undamped  members  to  find 


w 


(84a) 


E (t)  + (a+b\2)E  + A2(G*E  ) = R (t) 
m m m m m m 


where 


R (t) 
m 


- ,Ai  v > + < ,A2V  > 

^ L n 3 m s J 


(84b) 


and  where  the  proportionality  condition  in  equation  (30)  has  been 


employed.  Note  that  the  condition  B(’p  ,’■!>)  = 0 has  also  been  used. 

n m 


The  initial  condition  for  E (t)  may  be  found  by  premultiplying  the 

m 


initial  condition  form  of  equation  (73)  by  ^^2  and  then  integrating 


over  x to  give 


® E (0) 

< ibn, A u(x,0)  > = < ^n>A2Vs(x,0)  > + l < :^n, A2^m  > (85) 

m=l  m 


After  applying  undamped  orthogonality,  the  generalized  initial  conditions 
are 


ECO) 

d 


A2U 


(x)  > - < i)  ,A2V  (x 
u ms 


,0)  >J 


(86) 


and  similarly, 


Em(°) 

m 


VA2VS 


(x)  > - < * ,A2V  (x 
m s 


,0)  >J 


(87) 


An  alternative  form  of  the  right-hand  side  of  equation  (84a)  may 
be  found  which  incorporates  the  complex  compliance  J(t)  of  the  visco- 
elastic material.  With  this  form,  either  the  experimental  or  modeling 
information  of  a given  material  may  be  easily  included  in  the  solution 


for  E (t).  Begin  by  applying  the  proportionality  condition  in  equation 
m 


(30)  to  give 
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< ,AiV  > * (a+bA2)  < \p  ,Aov  > - b B(ip  ,v  ) 
m 1 s m m * s m s 


(88) 


Substitution  of  this  condition  into  equation  (84b)  leads  to 


1 [“• 


3+32)  < > A-2 vg  > - b B(iJ> 


wv] 


(89) 


where  R denotes  the  right-hand  side  of  equation  (84a),  and  c is  the 


m 


m 


proportionality  constant,  (a+A2b).  The  inner  product  mav  now  be  trans- 

in 

formed  by  use  of  the  Green's  identity  into 


< l ll  ,A2V  > » < v , Ay'll  > 
m * s s L m 


(90) 


Using  equation  (10a)  in  equation  (90)  gives 


< tb  ,A2 v > * - — < v ,Dib  > 
m ^ s x2  s’  m 

m 


(91) 


Now  apply  the  extended  Green's  identity, 


< v , Dip  >=<t(/,Dv  > + B(iji  ,v  ) 

s m ms  ms 


(92) 


to  obtain  from  equation  (91)  the  result, 


< ill  ,AyV  > m 
m s 


- h * *.'D%  > - h 

m m 


(93) 


The  term  Dvg  may  now  be  transformed  by  taking  the  Laplace  transform  of 
equation  (80) 


sGDv  * -F 
s 


(9/0 
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where  the  bar  represents  the  corresponding  transformed  variable  and 


G*f  = s G f (95) 

m 

Equation  (94)  may  be  easily  solved  for  Dvg  to  obtain 

Dv  = - — - -sJF  (96) 

S sG 

where  the  relation  between  the  Laplace  transforms  of  the  complex 
compliance  J and  the  relaxation  modulus  G has  been  used.  Taking  the 
inverse  Laplace  transform  of  equation  (96)  gives 


Dv  = - J * F 
s 


(97) 


Finally,  substitute  this  equation  back  into  equation  (93) , which  when 
placed  in  equation  (89)  leads  to 


R = - — 
m Qm 


f (cm3+32) 

< * , J * F > - BOP  ,v  ) 

^ X2 

in  ms 

v.  m 

• 

(98) 


-b  B(’p  ,v  ) 
m s 


This  is  the  alternate  form  for  the  right-hand  side  of  equation  (84a) 

which  was  sought.  This  form  of  R (t)  involves  directly  the  complex 

m 

compliance  J of  the  viscoelastic  material.  When  substituted  back  into 

equation  (84a)  the  solution  for  E (t)  may  now  be  found  by  taking  several 

in 

approaches.  If  a known  viscoelastic  model  composed  of  springs  and  dash- 
pots  has  been  selected,  the  compliance  will  be  known  once  the  model 
parameters  have  been  defined.  These  parameters  may  be  prescribed  a 
priori  or,  chosen  by  fitting  experimental  data  to  the  model  [14].  In 
addition,  experimental  data  may  be  used  directly  in  equation  (98)  without 
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r 


I 

I 

I 


any  reliance  on  a particular  model  [15],  in  which  case  equation  (84a) 

may  be  solved  by  Laplace  transforms,  which  may  be  inverted  analytically 

for  simple  models.  An  example  of  how  this  procedure  may  be  carried  out 

for  a Voigt-Kelvin  material  will  be  shown  subsequently.  Otherwise,  a 

numerical  inversion  process  may  be  needed.  Another  approach  which  may 

be  taken  in  finding  the  coordinates  E (t)  from  equation  (84a)  has  been 

m 

suggested  by  Valanis  [3].  When  a Laplace  transform  of  equation  (84a) 
has  been  taken,  a Voltera  integral  equation  of  the  second  kind  results. 
Then,  any  number  of  techniques  may  be  used  to  solve  this  integral  equa- 
tion. Note  that  if  the  usual  normal  mode  approach  had  been  undertaken, 

the  above  comments  on  how  to  find  E (t)  may  be  applied  unchanged  in 

m 

determining  C (t)  in  equation  (67).  However,  as  mentioned  previously 
the  acceleration  method's  more  rapidly  convergent  solution  makes  solving 
the  additional  quasi-static  problem  worth  the  trouble,  since  vg  may 
be  found  easily  for  many  structural  members  by  using  a correspondence 
principle. 

To  summarize  breifly,  the  dynamic  response  of  a viscoelastic 
structural  member  on  a proportional  viscous  foundation  may  be  formally 
expressed  by 

• E (t) 

u(x,t)  - v (x,t)  - l —r (x)  (99) 

s , u m 

m-1  mi 


where  ^(x)  and  are  the  corresponding  undamped  elastic  member  eigen- 
functions and  norm,  respectively.  The  temporal  coefficients  may  be 
found  by  solving 


E (t)  + (a+bX2)E  + \2(G*E  ) = R (t) 
m d m m m in 


(100) 
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and  R (c)  may  take  on  either  of  the  two  forms  as  given  by  equations 
m 

(84b),  or  (98).  As  special  cases,  when  a » b ■ 0,  a viscoelastic 
member  without  a damped  foundation  may  be  examined,  and  when  the  bilinear 
form  is  absent  the  member  possesses  homogeneous  boundary  and  in-span 
conditions. 

Forced  Response  of  a Structural  Member  with  Linear  Hvsteretic  Damping 

In  this  section,  the  concept  of  hysteretic  damping  will  be  explored 
for  the  purpose  of  finding  the  forced  response  of  structurally  damped 
structural  members.  The  hysteretic  damping  to  be  considered  is  linear; 
hence,  the  standard  assumed  solution  of  modal  analysis  will  again  be 
employed.  Both  self-adjoint  as  well  as  non-self-adjoint  systems  will 
again  be  considered. 

Consider  a continuous  structural  member  with  linear  hysteretic 
damping  governed  by  the  following  equation  of  motion: 

2 

Du(x,  t)  =»  £ A.  (x)3  .u(x,t)  - f(x)ei'‘t  (101a) 

j-1  3 3 

with  the  boundary  and  in-span  conditions, 

I^uU.t)  - P (x)ei0t  on  S (101b) 

where  f(x)  and  P^(x)  are  vectors  of  functions  which  represent  the 
spatial  variations  of  the  forcing  function,  and  boundary  and  in-span 
conditions,  respectively.  Note  that  the  temporal  portions  of  the 
excitation  functions  are  harmonically  oscillating  with  a driving 
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frequency  A,  which  is  a simplification  of  F(x,t)  given  by  equation  (la). 
The  initial  conditions  are  zero  because  interest  is  focused  on  the  steady 
state  response,  or  the  particular  solution  of  the  differential  equation. 

Based  upon  the  form  of  the  forcing  function  and  boundary  condi- 
tions, the  steady  state  response  sought  will  also  be  harmonic  with 
frequency  fl.  This  is  identical  to  a single  degree  of  freedom  system, 
since  all  points  in  the  member  are  oscillating  in  phase.  Hence 


(u,t)  - 3u  - iShi(x,t) 


(102) 


which  when  substituted  into  Equation  (101a)  yields 


Du(x,t)  ■ iflAju  + A2 3 2 u ~ f(x)ei“t  (103) 

This  matrix  equation  represents  the  governing  equations  of  motion  of  a 
structurally  damped  member  undergoing  steady  state  motion.  As  before, 
the  goal  is  to  determine  the  proper  modal  expansion  which  gives  the 
response.  This  entails  defining  the  modes  and  orthogonality  condition, 
and  the  uncoupling  coordinates.  It  will  be  shown  that  only  the  undamped 
modes  are  needed  to  resolve  the  solution  so  that  the  response  may  be 
written  as 


u(x,t) 


00 

l 

m-1 


q (x) 

m m 


(104) 


Hence,  the  coordinates  q (t)  must  be  determined  to  complete  the  steady 

IQ 

state  response. 
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Once  again  the  starting  point  of  the  general  formulation  is  the 


extended  Green's  identity, 

< ,Du  > - < u,D4»  > * B(ip  ,u) 

m mm 


(105) 


Since  Ax  and  A2  are  restricted  to  be  only  functions  of  x,  (or  at  best, 
spatial  matrix  differential  operators  with  homogeneous  boundary  condi- 
tions) B('|)  ,u)  contains  all  the  elastic  non-homogeneous  boundary  and 
m 

in-span  conditions  of  the  member.  Substitute  the  classical  free  vibra- 
tion problem  (equation  (3a))  and  equation  (103)  into  the  above  identity 
and  after  some  rearrangement  obtain 


^2  < ^m>A2u  < + iH  > iJ/^jAju  > + < iJ^,A2U  > 


< ii»  ,f(x)eiflt  > + B (ip  , u) 
m m 


(106) 


where  the  classical  form  of  equation  (25)  has  been  used.  Note  that  the 
imaginary  unit  i makes  the  operator  complex.  However,  by  foregoing  the 
Hermitian  norm  the  Green's  identity  remains  the  same  for  both  real  and 
complex  operators.  That  is,  the  imaginary  unit  is  effectively  treated 
as  a constant.  By  defining  ?m(t)  as  previously  (see  equation  (32),  the 
above  equation  may  be  rewritten  as 


l (t)  + X2?  + in  < <p  ,Axu  > 

in  mm  m 

■ < ,f(x)eint  > + B(J>  ,u) 
m m 


(107) 
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This  expression  is  not  yet  uncoupled  due  to  the  Ai  term.  A proportion- 
ality condition  must  be  invoked,  which  for  hysteretic  damping,  may  be 
called  proportional  structural  damping.  Recall  that  for  a single 
degree  of  freedom  system,  an  equivalent  viscous  damping  coefficient  can 
be  defined  for  hysteretic  damping  which  is  proportional  to  the  stiffness 
and  inversely  proportional  to  the  driving  frequency.  For  a continuous 
member,  a similar  equivalence  may  be  established  so  that 


Aju 


_ X 

Q 


Du 


(108) 


where  Meirovitch  refers  to  y as  the  structural  damping  factor  [16] 
Using  this  condition  leads  to 


if2  < ,Aiu  > = -iy  < i|>  ,Du  > = - iy  < u,D:f>  > 
• m m m 


-iy  B(i|)  ,u) 
m 


which  after  employing  classical  free  vibration  gives 


in  < ,Aiu  > * iX2y  < u,Ao<P  > -iyB(ty  u) 

m m m m 


iA2Y  < ’L’A2u  - iYB(*  ,u) 

mm  m 


(109) 


(110) 


Thus  equation  (107)  becomes 


k (t)  + X2(l  + iy)  5 - eint  < i Ji  ,f(x)> 

mm  m m 


+ (1  + iy)  B (*  ,u) 

m 


(111) 
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This  equation  is  identical  in  form  to  an  uncoupled  equation  for  a SDOF 
system  with  hysteretic  damping.  The  solution  for  non-homogeneous 
(time-dependent  or  constant)  boundary  and  in-span  conditions  may  be  found 
easily.  For  example,  if  the  boundary  and  in-span  conditions  are  separable 
such  that  the  right-hand  side  of  the  above  equation  may  be  written  as 


f eint  = [<  i>  ,f(x)  > + (1  + iy)  B (ij)  ,U(x)  ] eint 
mm  m 

f eint 

the  solution  for  5 (t)  is  given  by  = — 

" x2(i  + iy)  - n2 

m 


(112) 


(113) 


This  is  the  particular  solution  of  equation  (111)  since  it  is  assumed 
that  all  initial  condition  transients  have  died  out.  Since  classical 
orthogonality  was  used,  q (t)  in  equation  (104)  may  be  shown  to  be 


\(t) 


^m 


(114; 


so  that  the  steady  state  response  may  be  readily  found.  This  result  is 
identical  to  that  obtained  by  Meirovitch  [16],  except  that  this  formulation 
now  includes  non-self-adjoint  systems,  non-homogeneous  boundary  and  in- 
span conditions,  as  well  as  incorporating  the  very  powerful  matrix 
differential  operator  notation. 


Applications : 

In  this  section  will  be  presented  several  applications  of  the 
general  dynamic  theory  to  typical  damped  structural  members.  One  objec- 
tive of  presenting  these  examples  is  to  illustrate  in  detail  usage  of 
the  general  formulas  provided  in  previous  chapters.  However,  these 
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formulas  provide  the  dynamic  response  in  terms  of  some  quantities  which 

have  been  assumed  to  be  known.  For  example,  the  solution  of  the  damped 

free  vibration  problem  was  not  actually  solved  but  merely  represented 

symbolically  by  the  vector  of  eigenfunctions  v^  and  the  corresponding 

eigenvalue  S . Hence,  another  goal  of  these  examples  is  to  demonstrate 
n 

how  the  response  may  be  calculated  for  some  simple  damped  structural 
members.  It  will  be  shown  that  transfer  matrices  may  be  easily  coupled 
to  a modal  analysis  of  structural  members  which  possess  one  independent 
coordinate.  The  dynamic  response  is  computed  for  two  typical  members,  a 
beam  on  a non-proportional  viscous  foundation,  and  a viscoelastic  beam. 
Since  no  similar  results  could  be  found  in  the  literature,  these  cases 
will  be  presented  as* benchmark  examples.  For  more  complicated  members, 
the  response  will  be  theoretically  formulated  and  then  compared  with 
available  theoretical  results.  The  intention  here  is  to  show  how  the 
general  formulas  presented  in  this  section  lead  to  the  dynamic  response 
without  the  involved  algebraic  manipulations  taken  in  the  past.  Whether 
the  applications  are  numerical  or  theoretical,  these  examples  are 
intended  to  be  illustrative  rather  than  comprehensive.  Therefore,  com- 
plexities in  a structural  member  which  may  be  treated  without  confusion 
in  theory  but  are  tedious  in  practice,  will  be  avoided  to  demonstrate 
clearly  the  methodology  of  finding  the  dynamic  response. 

Cantilever  Beam  on  a Damped  Viscous  Foundation 

Consider  the  uniform  cantilever  beam  shown  in  Figure  1,  which  is 
resting  on  a stepped  damped  foundation.  The  value  of  the  mass  density 
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Figure  1 Damped  Cantilever  Beam 
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(p),  modulus  of  elasticity  (E) , moment  of  inertia  (I),  length  (£), 
and  viscous  foundation  coefficient  (c)  are 

p - .002301  lb  - sec2/in2 
E - 10  x 10s  psi 


I - .7854  inM 


l = 100  in 


(115) 


J.  0011505  lb-sec/in2  for  x < 50  in 

Y 0017257  lb-sec/in2  for  x > 50  in 

The  force  P applied  at  x=10Q  in.  is  a static  load  of  20  lbs.  At  t = 0+  it 
will  be  removed  and  the  previous  formulas  will  be  used  to  find  the  resulting 
transient  response. 

The  equations  of  motion  governing  this  member  have  been  introduced  in 
equations  (1).  The  matrix  differential  operators  are  given  by 


0 “fe 


1_  i_ 

3x  ~ El 


(116a) 


(116b) 
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(116c) 


- 


9 

M 


In  addition,  the  forcing  function  F(x,t)  is  set  equal  to  2ero  and  the 
initial  condition  vectors  are 


wo 

( 

o 

V 

uo(x)  = { 

o 

CD 

> V*’  ■ ( 

° 

0 

0 

0 

(117) 


where  w^  and  9q  are  the  static  deflection  and  rotation  due  to  P given  by 


w0  * 2EI 


(118) 


3 „ I*  ( * _ i) 

0 El  4 


(119) 


Note  that  Che  proportionality  condition  given  by  equation  (30)  is  not 

satisfied  by  the  variable  damped  foundation.  However,  this  stepped 

foundation  does  not  intefere  with  the  self-adjoint  property  of  this 

member.  Therefore,  the  damped  adjoint  eigenfunctions  v are  not  needed 

m 

to  determine  the  response.  For  this  case,  use  of  equations  (22),  (37), 
and  (45b)  leads  to  the  deflection  response. 


4. 

1 


w(x,t)  *^7  I 


2S  p ^ (w  )2dx  + (w  )2dx 
mm  m 


(120) 


ti. 

(c+S  p)(x2l  - x3/3)  dx  \w  (x) 
m / m 


In  this  expression  S is  the  complex  eigenvalue  and  w is  the  deflected 

m m 

mode  shape,  both  of  which  are  found  from  the  damped  free  vibration  prob- 
lem. This  free  motion  problem  may  not  be  solved  analytically.  Rather, 
a numerical  technique  must  be  employed.  Since  this  beam  possesses  one 
independent  coordinate  direction,  transfer  matrices  may  be  used.  For 
further  simplicity,  approximate  the  continuous  member  by  dividing  it 
into  massless  elastic  sections  punctuated  by  mass  lumps  (stations).  The 
damped  foundation  is  also  concentrated  into  discrete  dashpots  located  at 
the  mass  stations.  For  example,  this  leads  to  the  equation, 


(121) 


S^^+Sc^ 
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w 

9 

0 

“ [T  ] 

1 MJi 

M 

M 

V 

i 

V 

(122) 


where  the  superscripts  R and  L refer  to  locations  just  to  the  right  and  left, 
respectively,  of  the  mass  lump  (nO  and  dashpot  (c^)  at  x * x^.  The  transfer 
matrix  corresponding  to  the  mass  and  damping  lump  at  x = x^,  [ ] ^ , allows  the 
state  variables  just  to  the  right  of  the  mass  station  to  be  known  in  terms  of 


those  just  to  the  left.  In  a similar  wav,  the  expression  for  an  elastic 
section  is  given  by 


or. 
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(123) 


(124) 


where  A£^,  E^,  and  1^  are  the  length,  modulus  of  elasticity,  and  moment 

R L 

of  inertia,  respectively,  of  the  elastic  section  from  xi  to  xi+1  . 
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The  matrix  [T£]^  transfers  the  known  state  variables  at  x^^  to  state 
variables  at  x^+^.  With  equations  (122)  and  (124)  the  response  in 
terms  of  the  deflection,  slope,  moment,  and  shear  at  any  location  may  be 
found  in  terms  of  the  known  values  at  x ■ 0 by  merely  alternately  multi- 
plying mass  and  elastic  transfer  matrices  down  the  beam.  It  is  obvious 
that  variations  in  the  mass  density,  damping  foundation  modulus,  and 
moment  of  inertia  can  be  easily  accounted  for  by  merely  changing  m , 
c^,  E^,  and  1^  for  each  transfer  matrix.  Other  effects  such  as  rotary 
inertia,  shear  deformation,  in-span  conditions,  and  an  elastic  foundation 
may  also  be  readily  included  by  introducing  similar  transfer  matrices. 

Such  catalogues  of  transfer  matrices,  including  a more  thorough  treatment 
of  this  approach,  may  be  found  in  several  sources  such  as  Pestel  and 

Leckie  [17],  and  Pilkey  [18].  A listing  of  the  damped  frequencies  found 

« 

for  this  example  are  given  in  Table  1.  Using  the  damped  mode  shapes 
found  above,  the  response  may  be  computed  by  using  equation  (120).  This 
response  is  shown  in  Figure  2. 

In  conclusion,  a simple  Euler-Bemoulli  beam  has  been 
presented  to  demonstrate  how  easily  the  response  may  be  found  using  the 
general  formulas  presented  in  this  work.  The  response  may  be  found  for 
more  complicated  loading  conditions,  and  non-homogeneous  boundary  and 
in-span  conditions  by  using  the  same  equations.  The  approach  is  identical 
to  that  taken  for  this  simple  beam,  only  the  details  become  a bit  more 
tedious.  In  addition,  complications  such  as  shear  deformation  and 
rotary  inertia  are  also  readily  included  in  the  previous  analysis.  This 
is  done  by  modifying  the  appropriate  transfer  matrices  in  the  damped 
free  motion  problem  and  by  altering  the  matrices  A!  and  A2  used  in  the 
general  formulas . 
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Table  1 


Figure  2 Transverse  Tip  (x  = 100  in)  Deflection  of  Beam  on 
Variable  Viscous  Foundation 


Dynamic  Response  of  a Voigt-Kelvin  Beam 

Consider  the  undamped  Euler-Bemoulli  beam  shown  in  Figure  3,  which 
is  composed  of  a Voigt-Kelvin  viscoelastic  material.  Assume  that  the 
linearly  varying  distributed  load  shown  in  Figure  3 is  applied  at  t = 0 
and  then  maintained  at  that  value  for  all  future  time.  That  is, 

q(x,t)  = j H(t)  (125) 

where  H(t)  is  the  Heaviside  unit  function,  which  possesses  the  property 

H(t)  * J(t)  = J(t)  (126) 

The  dynamic  deflection  response  for  this  beam  is  sought  for  this  loading 
condition. 

The  governing  equations  of  notion  for  a general  viscoelastic  beam 
have  been  given  by  Robertson  [15].  The  equations  may  be  reduced  to 
one  fourth  order  differential  equation  when  shear  deformation  and  rotary 
inertia  are  neglected.  When  written  in  the  form  of  equations  (57) , the 
following  values  are  assigned: 


D - 2(1  » v)  I -2—  A,  - 0 

(127a) 

3x4 

A2  * -P  F » -q 

(127b) 

u - w(x, t)  uQ  - uQ  - 0 

(127c) 

'Vx-O.t  ' 1 a2)*.0,l  ’ 

(127d) 

■ <Vx-0,l  ■ 0 

(127e) 
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The  dynamic  response  of  this  viscoelastic  beam  may  be  found  from 
equation  (73)  to  be 


- E (t) 

w(x,t)  = w (x,t)  + l — ill  (x) 

s m*l  Tn  m 


(128) 


where  ws(x,t)  is  the  quasi-static  deflection.  For  this  simple  member, 
the  deflected  mode  shapes  'J'm(x)  and  corresponding  eigenvalues  X^  may 
be  analytically  determined  by  solving  the  undamped  eigenvalue  problem. 
That  is, 


mirx 


>li  (x)  * sin 

01  X* 


m * 0,±1,±2, 


(129) 


and 
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so  that  the  corresponding  norm  may  be  explicitly  determined  to  be 


Qm  ” P*/2  (131) 

Note  that  since  the  elastic  beam  is  self-adjoint,  only  the  modes  \p  (x) 

m 

are  needed  to  determine  the  response.  Therefore,  in  the  series  portion 

of  equation  (128),  only  the  temporal  coefficient  remains  to  be  determined. 

To  find  E (t),  recall  from  equations  (84a)  and  (98)  that 
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where  a * b = 0 since  the  beam  is  undamped.  In  addition,  the  initial 

values,  E (0)  and  E (0)  given  by  equations  (86)  and  (87),  are  zero 
m m 

because  the  initial  conditions  u^  and  u^  are  zero.  Introduce  equations 
(125),  (129),  (130),  and  (131)  into  the  above  expression  to  arrive  at 
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where  equation  (126)  has  been  utilized.  Eliminate  the  time  variable  in 
this  expression  by  taking  the  Laplace  transform  to  obtain 
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where  the  initial  conditions  of  the  complex  compliance  vanish  since  they 

are  proportional  to  E (0)  and  E (0)'.  Solve  for  the  Laplace  transform 

m m 

of  E (t)  (i.e.  E (s))  to  obtain 
m m 


-s3JJL 

E - 15 ' (135) 

m (A2  + s3J) 
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where  for  convenience 
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The  Laplace  transform  of  the  complex  compliance  for  a Voigt-Kelvin 
material  [19]  is  given  by 
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where 
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(138) 


* " qo/qi 


and  a.Q  and  are  Che  spring  and  dashpot  constants  in  the  Voigt-Kelvin 
model.  Introduce  equation  (137)  into  equation  (135),  and  after  some 
rearrangement  arrive  at  the  expression 
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where 
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and  where  the  method  of  partial  fractions  has  been  utilized.  Using  a 

table  of  Laplace  transforms,  E (t)  may  be  found  by  inverting  equation 

m 

(139)  which  leads  to 


E (t)  - -L 
m m 


X /X2q2-4q 
m mnl 


(rie 


r2t 


- r2e 


X 

m 


< 


/ 


Xmqi”4qc 
mi  c 


r2t 

(e 


-<t 

e 


] 


(141) 


This  is  the  expression  which  was  sought  for  E (t).  Since  it  can  be 

m 

evaluated  from  this  equation,  all  the  ingredients  needed  for  the  series 
portion  of  the  solution  in  equation  (128)  have  been  found  in  closed 
form.  To  determine  the  quasi-static  part  of  the  response  in  this  equa- 
tion, make  use  of  a so-called  correspondence  principle.  That  is,  as 
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Flugg  [20]  explains,  the  static  deflection  of  a viscoelastic  member  may 
be  found  by  finding  the  deflection  of  the  corresponding  elastic  member 
and  then  replacing  the  modulus  of  elasticity  (E)  by  the  ratio  of  poly- 
nomials Q(s)/P(s).  Here  Q(s)  and  P(s)  are  the  Laplace  transforms  of  the 
temporal  differential  operators  P and  Q,  introduced  in  equations  (135). 
In  general,  to  find  the  quasi-static  response  ws(x,t),  utilize  this 
correspondence  principle  for  each  time  step.  For  this  case  it  will  be 
shown  that  this  principle  need  only  be  applied  once,  since  the  load  is 
constant  for  all  time.  From  elementary  strength  of  materials,  the 
static  deflection  of  an  elastic  simply  supported  beam  is 


w(x)  - Y80"eT  [3x4  " 5x2£2  + 2*4] 


(142) 


Now  replace  q and  E by  (x/S, ) — and  Q(s)/P(s),  respectively,  to  obtain 
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which  is  the  Laplace  transform  of  the  quasi-static  deflection.  Note 
that  for  a Voigt-Kelvin  solid  [19], 
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equation  (143)  may  be  easily  inverted  to  give  for  the  quasi-static 
deflection  response. 
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which  decays  as  time  grows  since  k > 0.  Hence,  the  solution  of  a Voigt- 
Kelvin  beam  given  by  equation  (123)  has  been  determined.  It  is  apparent 
that  many  other  methods  of  solution  of  equation  (132)  could  have  been 
pursued.  For  example,  had  J(t)  been  experimentally  determined,  a 
numerical  solution  of  equation  (132)  would  have  been  attempted.  Or, 
this  experimental  data  could  be  fitted  to  a more  complicated  but  more 
realistic  model.  However,  for  more  complicated  viscoelastic  models,  a 
drawback  to  the  previous  method  of  solution  is  that  an  inverse  Laplace 
transform  of  a complicated  function  has  to  be  taken.  In  most  cases  of 
this  sort,  a numerical  inverse  Laplace  transform  would  probably  be 
warranted. 

Summary 

It  has  been  shown  in  this  work  that  a general  theory  for  the 
dynamic  response  of  linear  damped  continuous  structural  members  may  be 
formulated  using  a modal  analysis.  A general  set  of  formulas  was 
derived  to  provide  the  ingredients  needed  to  construct  the  dynamic 
response.  The  ingredients  included  a statement  of  the  free  vibration 
problem  and  a determination  of  the  orthogonality  relation  and  the 
uncoupled  temporal  coefficients.  As  shown  in  the  previous  chapters, 
these  elements  of  a modal  solution  were  found  in  terms  of  the  (assumed 
known)  governing  equations  of  the  member.  Therefore,  the  dynamic 
response  can  be  explicitly  written  out  by  merely  inspecting  the  govern- 
ing equations.  And  since  these  formulas  apply  to  self-adjoint  as  well 
as  non-self-adjoint  structural  members  with  or  without  non-homogeneous 
boundary  and  in-span  conditions,  the  dynamic  response  of  the  most  general 
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type  of  structural  member  may  be  found.  When  uncoupling  was  not 
possible,  another  advantage  of  this  general  approach  became  apparent  by 
exposing  the  mathematical  limitations  imposed  by  the  physical  model. 

This  work  was  divided  into  three  distinct  sections,  depending  on 

the  type  of  linear  damping  present  in  the  structural  member.  Viscous, 

viscoelastic,  and  hysteretic  damping  were  treated.  The  formulas  derived 

for  the  ingredients  in  a modal  expansion  are  summarized  in  Table  2.  All 

the  formulas  needed  for  each  form  of  damping  are  not  explicitly  shown  for 

the  sake  of  clarity.  The  lengthier  formulas  are  only  referenced  in  this 

table,  where  parentheses  denote  the  primary  equations  for  the  terms  in 

the  dynamic  response.  As  can  be  seen  from  this  table,  knowing  only  the 

solution  of  the  free  vibration  oroblem  (i.e.,  'i>  and  X or  v and  S ) 

’ n n n n 

the  dynamic  response  may  be  readily  determined.  These  formulas  are 
applicable  to  damped  (or  undamped)  self-adjoint  and  non-self-adjoint 
structural  members  with  arbitrary  loading,  non-homogeneous  boundary  and 
in-span  conditions,  and  initial  conditions.  In  all  these  cases,  the 
general  formulas  are  applied  in  exactly  the  same  way  to  find  the  dynamic 
response. 

In  summary,  a comprehensive  theory  has  been  presented  for  the 
dynamic  response  of  linear  continuous  damped  structural  members.  The 
formulation  depends  only  on  the  equations  of  motion  and  is  independent 
of  the  specific  structural  member.  Hence,  the  dynamic  response  of  a wide 
class  of  linear  structural  members  may  be  determined  by  the  same  theory. 
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Table  2 Formulas  for  the  Dynamic  Response 
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F5-120 

Office  of  Flight  Standards 
Tedere 1 Aviation  Agency 
Washington,  D.C.  2055) 

Chief,  Research  and  Development 
Maritime  Administration 
Washington,  D.C.  202 35 

Deputy  Chief,  Office  of  Ship  Coaecrt. 
Maritime  Adolnierrec ton 
Washington,  D.C.  20235 
Attn:  Mr  C.l.  Russo 

Atomic  Energy  Commission 
Div.  of  Reactor  Devel.  4 Toch. 
Garmantovn,  Maryland  20767 

Ship  Null  Reeeerch  Coanltteo 
National  Research  Council 
National  Academy  of  Sciences 
2101  Constitution  Avenue 
Washington , D.C.  20416 
Acta:  Mr.  A.R.  Lytle 
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UNIVERSITY  OF  VIRGINIA 
School  of  Engineering  and  Applied  Science 

The  University  of  Virginia's  School  of  Engineering  and  Applied  Science  has  an  undergraduate 
enrollment  of  approximately  1,000  students  with  a graduate  enrollment  of  350.  There  are  approximately 
120  faculty  members,  a majority  of  whom  conduct  research  in  addition  to  teaching. 

Research  is  an  integral  part  of  the  educational  program  and  interests  parallel  academic  specialties. 
These  range  from  the  classical  engineering  departments  of  Chemical,  Civil,  Electrical,  and  Mechanical  to 
departments  of  Biomedical  Engineering,  Engineering  Science  and  Systems,  Materials  Science,  Nuclear 
Engineering,  and  Applied  Mathematics  and  Computer  Science,  In  addition  to  these  departments,  there  are 
interdepartmental  groups  in  the  areas  of  Automatic  Controls  and  Applied  Mechanics.  All  departments  offer 
the  doctorate;  the  Biomedical  and  Materials  Science  Departments  grant  only  graduate  degrees. 

The  School  of  Engineering  and  Applied  Science  is  an  integral  part  of  the  University  (approximately 
1,400  full-time  faculty  with  a total  enrollment  of  about  14,000  full-time  students),  which  also  has 
professional  schools  of  Architecture,  Law,  Medicine,  Commerce,  and  Business  Administration.  In  addition, 
the  College  of  Arts  and  Sciences  houses  departments  of  Mathematics,  Physics,  Chemistry  and  others 
relevant  to  the  engineering  research  program.  This  University  community  provides  opportunities  for 
interdisciplinary  work  in  pursuit  of  the  basic  goals  of  education,  research,  and  public  service. 


